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RESUMEN 

El resumen sera traducido al espanol por los editores. 
We quantify the energy radiated through all the collisionally excited lines in a 
photoionized nebula which is permeated by temperature fluctuations. We assume 
that these correspond to hot spots which are the results of an unknown heating 
process distinct from the photoelectric heating. We consider all the effects of using 
a higher mean temperature (as compared to the equilibrium temperature) due to 
the fluctuations not only on each emission line but also on the ionization state of 
the gas. If this yet unknown process was to radiate a fixed amount of energy, we 
find that the fluctuations should correlate with metallicity Z when it exceeds 0.7 
solar. The excess energy radiated in the lines as a result of the fluctuations is found 
to scale proportionally to their amplitude t 2 . When referred to the total energy 
absorbed through photoionization, the excess energy is comparable in magnitude 
to t 2 . 

ABSTRACT 

We quantify the energy radiated through all the collisionally excited lines in 
a photoionized nebula which is permeated by temperature fluctuations. We assume 
that these correspond to hot spots which are the results of an unknown heating 
process distinct from the photoelectric heating. We consider all the effects of using 
a higher mean temperature (as compared to the equilibrium temperature) due to 
the fluctuations not only on each emission line but also on the ionization state of 
the gas. If this yet unknown process was to radiate a fixed amount of energy, we 
find that the fluctuations should correlate with metallicity Z when it exceeds 0.7 
solar. The excess energy radiated in the lines as a result of the fluctuations is found 
to scale proportionally to their amplitude t 2 . When referred to the total energy 
absorbed through photoionization, the excess energy is comparable in magnitude 
to t 2 . 
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1. INTRODUCTION 

The temperatures of photoionized nebulae are 
observed to be significantly lower when derived using 
recombination lines rather than from forbidden line 
ratios (e.g. Peimbert et al. 1995). This phenomenon 
has been ascribed to the existence of temperature 
fluctuations permeating the nebulae. Assuming that 
the fluctuations inferred by various authors (e.g. Pe- 
imbert et al. 1995, Esteban et al. 1998, Rola & 



Stasinska 1994) are caused by an additional albeit 
unknown heating agent (beside photoionization), we 
proceed to quantify the energy contribution which 
this unknown heating process must contribute to the 
total energy budget of the nebula in order to account 
for the much higher temperatures characterizing the 
collisionally excited lines as compared to those in- 
ferred from recombination lines (or nebular Balmcr 
continuum). In order to study the effects of arbitrary 
temperature fluctuations, we first describe the modi- 
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fications made to the multipurpose photoionization- 
shock code mappings ic (Ferruit et al. 1997). In 
Section H, we present photoionization calculations in 
which we consider different levels of fluctuation am- 
plitudes and quantify how they alter the global en- 
ergy budget of the nebula. A brief discussion is pre- 
sented in Section ||. 

2. THE ENERGY EXPENSE CAUSED BY 
NEBULAR HOT SPOTS 

We present a few definitions followed by the pro- 
cedure we adopt to implement the effect of temper- 
ature fluctuations in the code mappings ic. 

2.1. Definitions oft 2 and mean temperature Tq 
Following Peimbert (1967), we define the mean 
nebular temperature, Tq, as follows 



To 



Jv 



(1) 



in the case of an homogeneous metallicity nebula 
characterized by small temperature fluctuations; n e 
is the electronic density, T the electronic tempera- 
ture and V the the volume over which the integration 
is carried out. The rms amplitude t of the tempera- 
ture fluctuations is given by 



J v n 2 e (T-T ) 2 dV 
TZf v njdV 



(2) 



Note that we simplified the expression presented by 
Peimbert (1967) whose definition of t 2 differs in prin- 
ciple with each ionic species density ni while in the 
above equations we implicitely consider only ionized 
H (by setting n H + = n e ). Since t 2 in this paper is 
not an observed datum but an a priori global prop- 
erty of the nebular model, such differences are not 
important. 

The intensity of a recombination line is in general 
proportional to T a while for a collisionally excited 
line it is proportional to T^exp(—AE/kT) where 
AE is the energy separation of the two levels in- 
volved in the transition. In either case, p and a 
typically lie in the range —0.5 to — 1. 

2.2. Approximating fluctuations as hot spots above 

T 

Since the fluctuations' amplitudes t 2 derived from 
observations are much larger that that predicted by 
photoionization models (cf. Perez 1997), the solution 
to this inconsistency resides not in adding an addi- 
tional uniform heating/cooling term to the thermal 
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Figure 1. A numerical simulation of tempera- 
ture fluctuations consisting of hot spots character- 
ized by an amplitude t 2 = 0.005 (from Eq. §). The 
mean temperature To obtained from applying Eq. [I] 
is 9800 K (horizontal solid line). Also shown the rel- 
ative position of T rec and T^pm]43&3- The fluctua- 
tions' lower bound temperature (to be associated to 
T eq ) is 9500 K. 



balance equations (like heating by dust grain pho- 
toionization) since this would simly result in a uni- 
form raise of T and not in larger fluctuations. What 
is required to increase t 2 in models is that such pro- 
cess be operating in a non-uniform manner across 
the nebula and the picture which we propose is that 
of many hot spots created by this heating process. 

It is beyond the scope of this paper to consider 
any particular physical process to account for the 
fluctuations nor to model them in details. Our 
aim is limited to study the effect of ad hoc tem- 
perature fluctuations a la Peimbert on the energy 
budget under the assumption that these arise from 
this unknown extra heating mechanism acting non- 
uniformly across the nebula. 

In photoionization calculations, it is customary 
to define and use at every point in the nebula a 
local equilibrium temperature, T eqi which satisfies 
the condition that the cooling by radiative processes 
equals the heating due to the photoelectric effect. 
In our hot spots scheme, by construction, T eq cor- 
responds to the temperature floor above which take 
place all the fluctuations. For illustrative purposes, 
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we show in Fig. |l| a possible rendition of fluctuations 
(the solid wavy line) characterized by an amplitude 
t 2 = 0.005 (from Eq. [j]) and consisting of hot spots. 
Because T eq is a minimum in the distribution of T 
fluctuations, it defines the null energy expense when 
calculating the extra energy emitted as result of the 
hot spots. 

2.3. The determination of Tq 

The fluctuations are taken into account by 
mappings ic only in the statistical sense, that is by 
defining and using temperatures which are derived 
from Tq. Given a certain amplitude of fluctuations 
t 2 , we use everywhere a value for the mean tem- 
perature which is derived from the computed local 
equilibrium temperature (T eq ) and which takes into 
account that T eq is a lower extremum in the distri- 
bution of T fluctuations (consisting of hot spots). 
To define To, we found adequate using the following 
expression (based on the inverse of Eq. ||) 

T ^T e9 [l + 7(7-l)* 2 /2r 1/7 . (3) 

but where the optimum value of 7 has been in- 
ferred using numerical simulations of the hot spots. 
One simulation shown in Fig.[l] is characterized by 
T eq = 9500K_and t 2 = 0.005. By fitting 7 so that 
the value of Tq matches the numerically computed 
value of 9800 K, we obtain that 7 » —15. Within 
the deduced regime in which 7 <C — 1, the above 
equation varies slowly with 7. The purpose of this 
expresion is that it defines consistently To whatever 
the value of the equilibrium temperature calculated 
by MAPPINGS ic. Interestingly, 7 (or Tq) is invari- 
ant to simultaneously scaling up and down of all the 
hot spot amplitudes (which is equivalent to increas- 
ing or decreasing t 2 ). On the other hand 7 = — 15 
strickly characterizes a specific distribution of hot 
spot frequencies and widths. It is the correct value 
for fluctuations resembling those depicted in Fig. [j] 
but for other radically different distributions of hot 
spot widths and frequencies, 7 — —15 would only 
provide a first order albeit acceptable estimate of 
Tq. Despite this caveat, we estimate the uncertain- 
ties affecting the final determination of Theat. to be 
< 20%. 

2.4. Recombination processes 

For small fluctuations, the temperature can be 
expanded in a Taylor series about the mean Tq. In 
the case of recombination lines the intensity, I reci of 
a given line is affected by a factor 

7 rec ocf Q a [l + a(a-l)t 2 /2] . (4) 



This expression which can be used to compute indi- 
vidual recombination line intensities in the presence 
of small fluctuations, is equivalent to calculating the 
intensity (which is proportional to T a ) using instead 
the effective temperature T rec 

f rec = (r) 1 /" ~ T [l + q(q - l)t 2 /2]V« . (5) 

As shown by Peimbert (1995 and references 
therein), the temperature fluctuations have in gen- 
eral much less impact on recombination than on 
collisional processes which are usually governed by 
the exponential factor. For this reason, we adopt 
the simplification of considering a single value of 
a = —0.83 for all recombination processes (such a 
is the appropriate value for the H/3 line at 10000 K). 
This approximation will allow us to use a single tem- 
perature T rec when solving for the ionization balance 
of H, He and all ions of metals (equations in which 
enter recombination rates). 

If we consider the fluctuations drawn in Fig. [j] as 
example, we see that the mean recombination tem- 
perature Tree derived using Eq. || lies slightly below 
To. 

Interestingly, in the case of the ionization- 
bounded dustfrce models considered in this work, 
the adoption of one temperature or another when 
calculating recombination processes does not affect 
the global energy budget of the recombination lines. 
In effect, the total number of recombinations taking 
place across the whole nebula must equal the num- 
ber of ionizing photons produced by the UV source, 
whatever the value and behaviour of the tempera- 
ture. It is a self-regulating process: a hotter neb- 
ula [which results in slower recombinations rates and 
hence in a plasma containing less neutral H] of the 
type discussed below will simply turn out more mas- 
sive in ionized gas in order that the number of re- 
combinations remains equal to the same number of 
ionizing photons. The sum of all hydrogen recombi- 
nation lines intensities will remain unchanged. Also 
constant is the total amount of heat deposited in the 
nebula through photoionization. 

2.5. Collisional processes 
To compute the forbidden line intensities, we 
solve for the population of each excited state of 
all ions of interest assuming a system of 5 or 
more levels according to the ion. (In the case 
of intercombination, fine structure and resonance 
lines, we treat those as simple 2 level systems.) 
More specifically, when evaluating the excitation 
(oc T^exp[-AE i:j /kT]) and deexcitation (oc T^) 



4 



BINETTE & LURIDIANA 



rates of a given multi-level ion, each rate ij (pop- 
ulation) or ji (depopulation) is calculated using Tq 
(instead of T eq ) and then multiplied by the appro- 
priate correction factor, either 



This excess energy radiated in the form of collision- 
ally excited lines can be normalized respective to the 
total photoheating energy available. This defines the 
quantity T heat . 



in the case of excitation, or 

0^ = 1 + ^(^-1)^, (7) 

in the case of deexcitation. These factors result in 
general in an enhancement of the collisional rates in 
the presence of temperature inhomogeneities. They 
are adapted from the work of Peimbert et al. (1995) 
and were applied to all collisionally excited transi- 
tions. 

For the case of the particular rendition of the 
fluctuations depicted in Fig.Q, the mean collisional 
temperature characterizing the [Oiii]A4363 line (Pe- 
imbert et al. 1995) is shown by the horizontal dash 
line. Despite the large AEij involved in the emission 
of the line, the distance between this temperature 
and Tq is nevertheless smaller than that separating 
To from T eq . Therefore, in some instances the extra 
energy radiated through hot spots might depend as 
much on the distance separating To from T eq than 
on the above correction factors. 

Inspection of the line ratios calculated with 
mappings ic using Tb and the above correction fac- 
tors confirms (as expected) that the higher is AEV, 
the higher the line intensity enhancement (at con- 
stant i 2 ). It can be shown on the other hand that 
the far infrared lines (or any transition for which 
cxp[— AE^ /kT] sa 1) are less affected by the fluctu- 
ations (similarly to the recombination lines) and can 
even become weaker as a result of the fluctuations, 
in contrast to most collisionally excited lines. 

2.6. The energy radiated through hot spots 
Our aim is to quantify the excess energy gener- 
ated by temperature fluctuations under the assump- 
tion that these are caused by a putative heating 
mechanism which operates within small regions ran- 
domly distributed across the nebula. To calculate 
this energy we simply integrate over the nebular vol- 
ume V the luminosity of each line (or transition) ij 
using the statistically determined local To and the 
multiplicative correction factors of Eqs || and 0, and 
then subtract the corresponding luminosity obtained 
by using the equilibrium temperature T eq instead. 



r = ZiiU^jjj-^ffl dv 

J2 l3 f v ^dV + E k J v €dV 

(8) 

Lfiuc — L eq 
~ T +0 ' ( ' 

where j{j U ° corresponds to the local nebular emis- 
sivity of line ij calculated using To and taking into 
account the above correction factors while j^ is the 
corresponding emissivity assuming equilibrium tem- 
perature everywhere. The term with q^ q corresponds 
to various cooling rates from processes not involv- 
ing line emission such as freefree emission, while Q eq 
represents the total volume integrated value of this 
term. L fi uc and L eq correspond to the integrated en- 
ergy loss across the whole nebula due to collisionally 
excited lines with fluctuations and without, respec- 
tively. Since T eq must satisfy the condition that the 
cooling rate equals everywhere the heating rate, we 
obtain that L eq + Q eq is also equal to the total energy 
deposited into the electronic gas by the photoelectric 
effect. 

A larger fraction of the ionizing radiation simply 
keeps the nebula ionized (resulting in recombination 
lines) but does not affect the nebular temperature. 
An alternative way therefore of expressing the im- 
portance of the excess cooling due to the fluctuations 
is to use as reference the total energy absorbed by 
the nebula including the energy emitted as recombi- 
nation lines and nebular recombination continuum. 
We define T a i, s . as follows 

r Lfiuc — L eq 

l abs. — yji , (1UJ 

^total 

with Lf otal the ionizing luminosity of the exciting 
star 

/>OC 

Lf otal = / L v dv , (11) 

Jvq 

where L„ is the energy luminosity distribution of the 
ionizing star and vq the frequency corresponding to 
the ionization treshold of H. To be consistent with 
this definition, we must consider only nebular models 
which arc ionization-bounded and fully covering the 
ionizing source (over 4-7T sterad). Depending on the 
hardness of the ionizing radiation, r a fc s . turns out 
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to be 2-3 times smaller than Theat. because of the 
larger fraction of the absorbed energy which goes 
into photoionizing rather than into heating the gas. 

In summary, the modifications made to 
MAPPINGS ic to consider the effects of T fluc- 
tuations not only include the calculation of the line 
intensities using the formalism described in Pcim- 
bert et al (1995) but also considers their impact on 
the ionization balance across the nebula through the 
use of T rec (instead of T eq ) to derive the recombina- 
tion rates. Since the nebula is substantially hotter 
on average when t 2 > 0, it will be more ionized 
since the recombination rates are slower at higher 
temperature. This in turn results in a lower pho- 
toionization and heating rates given that there is less 
neutral H in the nebula, therefore the equilibrium 
temperature T eq computed by MAPPINGS ic will be 
lower than without fluctuations. All these effects 
have been taken into account self consistently and 
do not affect the energy conservation principle in 
the case of ionization-bounded nebulae as discussed 
above and in in Sect. 
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3. MODEL CALCULATIONS 

We have explored the behaviour of Theat, in pho- 
toionization models of different metallicity (Z), ex- 
citation (U) and different spectral energy distribu- 
tions (hereafter SED). We will express the nebu- 
lar metallicity with respect to the solar abundances 
(from Anders & Grevesse 1989) for which we take 
that Z = 1. To define other metallicities, we sim- 
ply scale the abundances of all the metals respective 
to H by a constant multiplicative factor equal to Z. 
We summarize below our results under various model 
conditions. 

For the photoionized Hn regions, we have se- 
lected unblanketed LTE atmosphere models from 
Hummer & Mihalas (1970) of temperatures T e // of 
40000 K, 45000 K and 50000 K [see Evans (1991) for a 
comparative study of nebular models using different 
model atmospheres]. To represent planetary nebu- 
lae, we simply employed black bodies of 10 5 K and 
10 5 3 K truncated at 54.4 eV. The geometry adopted 
in the calculation is plane-parallel with a gas density 
of uh = 10 cm~ 3 in all cases. The excitation of the 
nebula is defined by the excitation parameter U as 
follows 



U 



1 



cn H 



Anr 2 hv 



dv = 



cn H 



(12) 



where c is the speed of light, h the Planck constant 
and r the distance of the slab from the ionizing star. 
U is the ratio between the density of ionizing photons 
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Figure 2. Behaviour of Theat. in a sequence of 
photoionization models which have different nebu- 
lar metallicities (relative to solar Z = 1). Two SED 
used in the calculations consisted of a 45000 K and a 
100000 K star, respectively In all models t 2 = 0.04 . 



impinging on the slab {lph/c) and the total H den- 
sity. All the calculations carried out were ionization- 
bounded. 

Since hotter SEDs result on the other hand in 
much higher photoelectron energies and hence hotter 
nebulae at all metallicities, the maximum of the long- 
dashed curve corresponding to the 10 5 K blackbody 
in Fig.|| is shifted towards higher Z relative to the 
45000 K SED. 

3.1. Dependence of excess heating on metallicity 
We have calculated nebular models of different 
metallicities covering the range 1% solar (Z = 0.01) 
to 4.7 times solar. In the models presented in Fig. || 
all other parameters are identical, namely, U — 0.01, 
t 2 = 0.04 and a spectral energy distribution (SED) 
having either T eff = 45000 K or 100000 K. It can 
be seen that a maximum in Theat. occurs within 
the range Z ~ 0.2-0.4. For the Hn region mod- 
els, the average values for T eq across the nebulae 
are ~ 15000 K, 9000 K, 5000 K and 1500 K for the 
Z = 0.01, 0.7, 2.5 and 4.7 models, respectively. 
Given that the nebular temperatures decrease mono- 
tonically with increasing Z, the curves' behaviour 
can be understood as follows: at very low metal- 
licities Z <C 0.2, the forbidden lines of metals are 
not the main cooling agent and the fluctuations have 



G 



BINETTE & LURIDIANA 



therefore a negligible impact on the total cooling. At 
higher Z values around solar, however, the cooling 
due to metals is very large and the nebula turns out 
much cooler, to the extent that many optical lines 
become now less intense despite the increase of the 
metal abundances. At even higher Z, the optical 
lines are progressively 'switched off' and cannot con- 
tribute to the cooling of the nebula. In this regime, 
these lines have AE <§C kTo and become somewhat 
brighter as the temperature is further lowered, ex- 
plaining why the fluctuations now cause Theat. to 
become negative. The final rise at the upper Z end 
in the H II region model sequence (solid line) reflects 
the fact that at such low temperature the cooling can 
only be carried out by the low energy transitions (far 
infrared lines) where again AE > kTo- 

The abrupt decrease of Theat. above solar metal- 
licity has interesting consequences if the unknown 
heating process responsible for the fluctuations 
needed to radiate a comparable amount of energy 
in different objects. This could arguably be the case 
for instance if this process was reconnection of mag- 
netic field lines. All other factors being equal, nebu- 
lae three times solar would require much larger am- 
plitude turbulences (larger t 2 ) to radiate the same 
amount of energy than a solar metallicity nebula. In 
this case we might expect to observe a correlation 
between t 2 and metallicity (beyond Z > 0.7). 

3.2. Dependence of excess heating on t 2 

In the following calculations, we adopt two repre- 
sentative metallicities of Z = 0.2 and Z = 1 (solar). 
We found no clear trends across different SEDs of 
how Theat. varied with U and therefore we only re- 
port results concerning a single ionization parameter 
of value 10 -2 . Campbell (1988) has shown that the 
range of U for most Hll galaxies lies in the range 
10" 2 - 6 to lCT 1 - 8 . 

In Figs.|| and [|, we show the behavior of both 
Theat. and r a fc s . as a function of increasing t 2 of the 
models. Each line corresponds to a given SED and 
Z. The increase in Theat. is steeper at small values 
of t 2 <J 0.1 followed by a more linear regime with a 
slope < 2.5 at large t 2 . The radiated energy con- 
tribution from the fluctuations (compared to photo- 
electric heating) is of order 2.5£ 2 (Fig. |J) but with a 
wide dispersion when Z = 1 (curves labelled S). For 
Z = 0.2 there is no dependence on the SED and the 
curves almost superimpose each other. 

If we consider the total energy budget (photo- 
heating plus recombination energy, see Eq. |l0|) and 
not just the photoelectric heating part, we find that 
Tabs, is comparable in magnitude to t 2 (Fig. ||) but 
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t 2 

Figure 3. Behaviour of Theat. with increasing t 2 . 
Two metallicities were used: solar (labelled S) and 
0.2 solar (labelled L) and 4 SED: 40000 K,45000K, 
50000 K, 100000 K and 200000 K labelled 40, 45, 50, 
100 and 200, respectively (the L sequences are too 
close to allow labelling of each SED). Theat. becomes 
approximately linear at larger t 2 . For comparison we 
show a thick line representing Theat. — 2.5i 2 . 
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t 2 

Figure 4. Behaviour of r a fc s . with increasing t 2 . 
Same notation as in Fig. ^. For comparison we show 
thick lines representing T aos , = t 2 and 2.5i 2 , respec- 
tively. 
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again with a wide dispersion which result from dif- 
ferences in either Z or the SED. 

4. DISCUSSION 

If the turbulences are the result of heating by a 
yet undiscovered process, we infer from Fig. |^ and ^| 
that the extra energy radiated via all the collisionally 
excited lines due to the turbulences is a substantial 
fraction (r a fc s . > 10%) of the total photoionization 
energy whenever t 2 J> 0.02 for planetary nebulae 
and t 2 ;> 0.08 for H n regions, respectively. Taking 
instead as reference only the heating by photoioniza- 
tion (T a b s ,), these limits are reduced by two. In this 
case almost all nebulae would radiate more than 10% 
(r o f, s .) of their energy as a result of T fluctuations. 
In effect, typical empirically determined t 2 values for 
galactic H n regions lie in the range 0.02-0.06 (Lurid- 
iana 1999). For extragalactic Hn regions (Luridiana 
et al. 1999, Gonzalez-Delgado et al. 1994) and plan- 
etary nebulae (Peimbert et al. 1995), even larger 
values have been encountered, of order 0.1 or more^]. 
In those cases, the energy envolved can be a large 
fraction of the energy budget especially if the excit- 
ing stellar temperature exceeds 10 5 K. 

The underlying assumption behind our calcula- 
tions is that an external heating agent is at work to 
generate T fluctuations. We should point out that 
there exist alternative explanations to the fluctua- 
tions which rest on photoionization alone and have 
not been completely ruled out. Possible mechanisms 
which we plan to study in some details: 

1. Metallicity inhomogeneities Temperature fluc- 
tuations would result naturally from nebu- 
lae consisting of a multitude of condensations 
of greatly varying metallicities (see Torres- 
Peimbert, Peimbert & Peha 1990) 

2. Gas expansion Fluctuations could arise from 
an outflowing wind generated at the surface of 
dense condensations (problyds) as a result of the 
champagne effect. Rapid adiabatic expansion 
would result in overcooled emission regions in 
the wake of the wind (resulting in negative T 
fluctuations with respect to T eq ). 

3. Ionizing source variability A rapidly varying 
ionizing field will generate a sequence of out- 
wardly propagating ionization and recombina- 
tion fronts. The basic asymmetry existing be- 
tween ionization fronts (propagating at a speed 

1 Measured values of t 2 in excess of 0.1 should be considered 
only as rough estimates since the assumed regime of small 
fluctuations does not apply anymore. 



limited by changes in opacity) and recombi- 
nation fronts (non-propagating), would lead to 
temperature fluctuations. 

To our knowledge, none of these mechanisms 
have yet been incorporated explicitly in any nebular 
model and we therefore cannot assess their particu- 
lar merits. Mechanisms which are good candidates to 
explain the extra heating assumed in this work (see 
also discussion by Stasihska 1998) and which are not 
directly related to photoionization, have been pro- 
posed and include the following processes: 

1. Shock heating Either as a result of stellar wind 
or from supernovae. In the case of giant H n 
regions, this later mechanism was favored by 
Luridiana et al. (1999). 

2. Magnetic heating Turbulent dissipation from 
Alfen waves and magnetic line reconnection has 
been proposed to explain the variations of the 
[N n] /Ha ratio emitted by the warm diffuse ion- 
ized gas in the Galaxy (Reynolds et al. 2000). 
In the case of planetaries, the large values of 
t 2 found would then make sense in the light of 
the success of the magnetically accelerated wind 
models of Garci'a-Segura (2000). Furthermore, 
Peimbert et al. (1995) has shown that the nebu- 
lae with the highest gas velocity dispersion show 
the highest values of t 2 , which would be consis- 
tent with an increasing role played by magnetic 
acceleration of the gas in the objects with the 
largest velocity dispersion. 

In summary, the current work shows how T a b s . 
and Theat. vary with metallicity and t 2 as well as how 
these quantities are affected by the ionizing energy 
distribution. Our results can be translated into en- 
ergy requirements (function of t 2 ) which competing 
explanations of the temperature fluctuations must 
satisfy and, hence, can be used to probe their re- 
spective viability. 

The work of LB was supported by the CONACyT 
grant 32139-E. 
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